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Abstract. Phonons in expanding Bose-Einstein condensates with wavelengths much 
larger than the healing length behave in the same way as quantum fields within a 
universe undergoing an accelerated expansion. This analogy facilitates the application 
of many tools and concepts known from general relativity (such as horizons) and the 
prediction of the corresponding effects such as the freezing of modes after horizon 
crossing and the associated amplification of quantum fluctuations. Basically the same 
amplification mechanism is (according to our standard model of cosmology) supposed 
to be responsible for the generation of the initial inhomogeneities - and hence the 
seeds for the formation of structures such as our galaxy - during cosmic inflation (i.e., 
a very early epoch in the evolution of our universe). After a general discussion of the 
analogy (analogue cosmology), we calculate the frozen and amplified density-density 
fluctuations for quasi-two dimensional (Q2D) and three dimensional (3D) condensates 
which undergo a free expansion after switching off the (longitudinal) trap. 
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1. Introduction 

A quarter of a century ago, W. G. Unruh P noticed an intriguing analogy between (at 
a first glance) very different systems: Phonons in irrotational fluids behave in the same 
way as (quantum) fields in curved space-times whose geometry is determined by the 
effective metric and depends on the velocity of the fluid flow and the speed of sound etc. 
This analogy facilitates a know-how transfer in both directions: On the one hand, it 
allows us to apply all the concepts, tools, and effects known from general relativity - such 
as horizons - to the propagation of phonons (or other quasi-particles) in fluids and hence 
leads to a better understanding of condensed-matter phenomena in terms of universal 
geometrical concepts. On the other hand, this analogy opens up the opportunity of 
a theoretical and possibly experimental investigation of exotic quantum effects known 
from cosmology such as Hawking radiation 

Unfortunately, it turns out that an experimental verification of the Hawking effect 
by means of such an analogue system (i.e., a black hole analogue) is very difficult 
jsJj. In this article, we shall focus on another exotic quantum phenomenon known 
from cosmology which also involves the amplification of the initial quantum vacuum 
fluctuations due to the presence of a horizon: Within our standard model of cosmology, 
the early universe was nearly homogeneous and basically all inhomogeneities - including 
the seeds for the formation of large-scales structures such as our galaxy - originate 
from the quantum vacuum fluctuations of a scalar field, which is called the inflaton. 
During inflation ([I]), which is a (conjectured) period of accelerated expansion of the 
very early universe, these quantum fluctuations were amplified due to the presence of a 
cosmic horizon (cf. Section EJ). These amplified fluctuations left an imprint in the cosmic 
microwave background radiation - small perturbations upon a homogeneous background 
as measured by the WMAP satellite jSj. 

Owing to the aforementioned analogy, the very same amplification mechanism 
should also occur in appropriate fluids. As an example, we shall consider the in some 
sense best understood superfluids - dilute atomic/molecular Bose-Einstein condensates 
(BECs), see, e.g., for review. Apart from the good theoretical understanding, further 
advantages of BECs are the various possibilities to manipulate (e.g., changing the speed 
of sound) and control them experimentally. As it will become evident later, the free 
expansion of a BEC after switching off the trap - which is a standard procedure in 
time-of-flight experiments/measurements, for example - leads to the formation of an 
effective cosmic horizon and, consequently, to the amplification of the initial quantum 
fluctuations of the phonon modes. These amplified fluctuations manifest themselves in 
potentially measurable small-scale density (and phase) variations. 

In order to study this effect for rather general conditions, we shall consider 
BECs with non-specified power-law self-interactions in an arbitrary number of (spatial) 
dimensions. After a discussion of the effectively lower-dimensional behaviour of strongly 
constrained BECs in Section|2l the effective metric governing the propagation of phonons 
is derived in Section |3] for arbitrary power-law self-interactions in any dimension. 
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Section |U is devoted to the introduction of co-moving coordinates (in complete analogy 
to cosmology) and to the scaling behaviour of the expanding condensate. The analogy 
to cosmology is further elaborated in Section El by applying the concept of an effective 
horizon, whose existence generates the aforementioned amplification mechanism. The 
spectrum and magnitude of the resulting density fluctuations are calculated in SectionsIB] 
and [7| for the quasi-two and the three-dimensional case, respectively, assuming the usual 
quartic coupling. 



2. Dimensional Reduction 



Since we shall consider an arbitrary number of spatial dimensions later on, let us 
first discuss the behaviour of strongly confined and hence effectively lower-dimensional 
Bose-Einstein condensates. In three spatial dimensions, Bose-Einstein condensates are 
described by the Lagrange density (h — 1 throughout) [Zj 

C = _ ***) - i^-L - V cxt (r,m\ 2 - I |v&| 4 , (1) 

with m being the mass of the bosons, V ext the external and generally space-time 
dependent one-particle trapping potential, and the two-particle coupling g in s-wave 
approximation, which is related to the s-wave scattering length a s via g = ATra s /m. If 
the external one-particle trapping potential can be split up into a parallel and a static 
transversal part via 

V c Ar,t) = vl t (r h t)+V e i t (r ± ), (2) 

it is useful to decompose the order parameter \1> (and its quantum fluctuations 5^/) into 
a complete set of real and time-independent functions a (r_|_) governing the transversal 
dependence (cf. where the density of an elongated BEC is decomposed in a similar 
manner) 

tt(r,t) = $>a(rj^ a (r,|,t). (3) 

a 

The time evolution of the coefficient functions if> a (r\\,t) is determined by the reduced 
Lagrangian density 

£n / ' dV ? r D C 

' ' 3-D 



r3-D 




= E - 4m - • - vlcvJ / dv_ 

a/3 J 

- | E / dVt D tidrftyi . (4) 

a(3-yS 

If we now choose the functions <p a {r_[_) to be the orthonormal eigenfunctions of the 
self-adjoint operator /C with eigenvalues Q a 

kk = ( _ S + v *) k = na<pa ' (5) 
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the reduced Lagrangian density simplifies to 

A = E 



i 1 



_ 2 TQT " / 2m 1 

- § E ^im>> / ^^l 3_d &<m> • (6) 

a[3-yS 

The last term induces a coupling of different modes, which complicates the solution. 
However, if we assume that the coupling term is sufficiently small and the lowest mode 
a = dominates \ipo\ 3> |?/Wo|, we may estimate the population of the higher modes in 
analogy to stationary perturbation theory: The mixing between the lowest mode a = 
and all higher modes a > is small if the energy differences AQ a = Q a — Qq are large 
compared to the transition matrix elements of the interaction Hamiltonian, i.e., 

AQ a » (# | 2 | dVt D <f>* a <f>f3<f>;<f>6- (7) 

Since the modes <fi a are normalised <p a ~ 1/ \Jvf_~ D , this condition can be re-expressed in 
terms of the size of the transversal dimension a± and the healing length £ = 1/ \J g\ipQ\ 2 m. 
Therefore, imposing the conditions 

£ > a± > a s , (8) 

where the latter requirement a±_ ^> a s is necessary for the Gross-Pitaevskii Lagrangian 
density in Eq. (JTJ) to be valid, we arrive at the reduced Lagrangian density for the lowest 
mode ip = ip 

A = \{Pi> - W) - W| 2 - vlM - 1 H 4 , (9) 



gi=g I dVt D |0 O | 4 = — / dVt D |0 O | 4 oc , (10) 



where Qo has been absorbed into V^t and 

4ira 

m J - a '^ 

denotes the reduced (lower-dimensional) coupling constant. For example, the 
eigenfunctions 4> a (f±) = <fia(z) of a harmonic potential V±(z) = mu 2 z 2 /2 with the 
transversal length scale a± = 1/ y/muj z are just the Legendre polynomials with an 
equidistant energy spectrum. In this case, the reduced coupling is given by 



#U = gymLU z /2n. However, it should be emphasised that the method presented above is 
applicable to more general potentials satisfying the aforementioned conditions as well. 

Note that if we relaxed condition (JEJ) and the transverse size of the condensate, a z , 
would be of comparable to the scattering length, a s /a z = 0(1), we would obtain a lower- 
dimensional condensate with different interactions due to density-dependent corrections 
to the coupling coefficient EH QU Q2! ■ 

3. Effective Geometry 

After having discussed the dimensional reduction, let us start with the Lagrange density 
in an arbitrary number of spatial dimensions D (where we omit the superscript || for 
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the sake of conciseness). In addition, we shall assume a more general self-coupling term 
\ip\ 2N which will be motivated later 

c = i^v - r^p) - ^1 wi 2 - v; xt H 2 - 1 \^\ 2N . (ii) 

Inserting the Madelung representation 

4> = ^e iS , (12) 
the Lagrange density reads 

C = -Qd t S - ^{VSf - - V ext g- 9 -g N . (13) 

2m 2m 2 

As usual, variation w.r.t. S yields the equation of continuity and w.r.t. g the Bernoulli 
equation with the mean-field velocity v = V^/m and the specific pressure p(g, V 2 ^). 
Assuming that the density profile is sufficiently smooth, i.e., that the relevant length 
scales are much larger than the healing length (Thomas- Fermi approximation), we 
neglect the quantum pressure term (V^) 2 . Linearisation S = S + SS and g = g + 5g 
around a background solution Sq and £o yields the second-order Lagrange density 

(cf. omnia eh) 

£( 2 ) = -6 e d t 6S - —CV5S) 2 -6gv - V6S -—6q*, (14) 
2m 2 

where we have introduced the effective coupling 

N(N-l) N _ 2 

9N = g g • ( 15 ) 

For the usual case = 2, we have = g, but one should still bear in mind that we 
are considering an arbitrary number of spatial dimensions D. Note that there are no 
sound waves at all for iV = 1, since the theory is non- interacting in that case. 

In the Thomas- Fermi approximation, the linearised Bernoulli equation can be solved 
for the density fluctuations 

5g = - 9t + V °- V SS, (16) 

and inserting this result back into Eq. (|14j) . we obtain the effective Lagrangian for the 
phase fluctuations <fi = 5S only 

4? = _i_(^ .V^-|l W )S. (17) 

Even in an arbitrary number of spatial dimensions D (the difficulties for D = 1 will 
be discussed below) and for general self-coupling N > 1, this Lagrangian is completely 
equivalent to that of a free (minimally coupled) scalar field in a curved space-time (e.g., 

P2I) 

l 



45 = 2vW(^)«S(^)> ( 18 ) 

provided that we insert the Painleve-Gullstrand-Lemaitre (PGL) metric jTHj 

^cfr _ AN) ( c%- vl v \ 
g ^"^ D { v -1 ) ' 

gZ = ^pJ N {v v ®v °-c%l)' (19) 
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with the speed of sound c 2 N = gNQo/m and the conformal factor 

, \ 2/(0-1) 

= p) • (20) 

\9nJ 

This expression already indicates problems in one spatial dimension D = 1 due to 
conformal invariance of the scalar field action (JTSJl in 1+1 dimensions. Without the 
introduction of an additional dilaton field, the identification of an effective metric is 
only possible if c^/gjy = const, for example if g = const and N = 3, see also the next 
Section. 

4. Co-Moving Coordinates and Scaling 

Since phase fluctuations in Bose-Einstein condensates behave (in the Thomas-Fermi 
approximation) exactly as a scalar field in a specific curved space-time described by the 
effective metric in Eq. (119)) . it will be useful to investigate this metric with the tool 
known from general relativity. The effective line element reads 

ds 2 cS = ([c 2 N - v 2 Q }dt 2 + 2v ■ drdt - dr 2 ) . (21) 

In analogy to cosmology, we shall assume local isotropy and homogeneity - which will be 
a good approximation in the centre of the BEC cloud. Since the coupling g is supposed to 
be constant for simplicity, this assumption implies a spatially homogeneous but possibly 
time-dependent density g = g (t) and effective coupling g N = gi<[{t). Furthermore, after 
a suitable re-definition of the origin of our coordinate system, we may set Vq oc r due to 
the presumed local isotropy and homogeneity. Insertion of this ansatz into the equation 
of continuity yields ([T4*| ITHj) 

^ o(t) = ^ir ~ v ^ r) = \^ (22) 

which allows us to describe Qo(t) and v (t,r) as well as #Ar(t) by means of a single 
time-dependent scaling parameter b(t) satisfying the initial condition b(t = 0) = 1 and 
b(t = 0) = 0. As we know from general relativity, an off-diagonal metric such as in 
Eq. ()21j) can be diagonalised by introducing co-moving spatial coordinates via (cf. the 
scaling transformation for BECs in [TTH I2U)) 

dp = dr ~ V ° dt - ds 2 eS = Ag° (c 2 N dt 2 - b 2 dp 2 ) . (23) 



r m m 

In addition, we may transform from the laboratory time t to the effective proper (co- 
moving wrist-watch) time r in order to eliminate the factor in front of dr 2 



r = jdt yjA^P c N -> ds 2 cS = dr 2 - A ( ^b 2 dp 2 , (24) 

arriving at the standard Friedmann- Robert son- Walker representation, see, e.g., ^7]. Let 
us investigate the remaining factor A^b 2 a bit further: Since Qo{t) oc b~ D (t) according 
to Eq. (J22J), we obtain g N (t) oc b- D{ - N ~ 2 \t) from Eq. JTSJ as well as c N (t) oc fe-- ^ -1 ^ 2 ^), 
i.e., the factor in front of dp 2 scales as 

A^b 2 OC b ^D(N-3)/(D-l) ; (25) 
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0, (27) 



^=(-i^ + Kxt(*, r) + g ($t)"-i$"-i $ ; (28) 



according to Eq. (|2(J|) . Interestingly, the exponent vanishes for 

N =J} + 1 > (26) 
e.g., for D = 2 and iV = 2, or D = 3 and iV = 5/3, or D = 1 and AT = 3. As we have 
observed above, special care is required for the latter case D — 1, but for N = 3 we 
can indeed introduce an effective metric and the conformal factor can be chosen at will 
since it does not enter the calculation. 

If the condition in Eq. ()26|) is satisfied and thus the exponent in Eq. ()25|) vanishes, 
the transformed metric in Eq. (|2"2|) is flat and hence the wave equation for the phonon 
modes becomes trivial in terms of the coordinates r and p 
f d 2 1 d 2 s 

i.e., a solution <p(r, p) is independent of the external time dependence mediated via the 
scaling parameter b(t). It turns out that this perfect scaling is not only valid for the 
phonon modes, but can be extended to the full field operator: If we start from the 
equation of motion in the Heisenberg picture 

d * fid 2 

, — \l> = „ 

dt \ 2m dr" 1 

again with general s-wave coupling (AT) and in D spatial dimensions, we may account 
for an arbitrarily time- dependent external one-particle trapping potential V ext {t,r) as 
long as it is purely harmonic at all times 

V ext (t,r) = -u; 2 xt (t)r 2 , (29) 

by inserting the scaling ansatz derived above (cf., [2*U] ) 
3 exp^rm^ r)/2} ^ 

The numerator expjzmv^t, r)/2} of the pre-factor reproduces Vq and the denominator 
^Jb D (t) accounts for Qo(t) = Qo(t = 0)/b D (t) and ensures the correct commutation 
relations for ^(r, p), i.e., 

mr,p),^(r,p')] = 5 D (p-p'). (31) 
For a flat metric N = 1 + 2/ D, the proper time is (independently of D) determined as 

wr (32) 

and p is of course still given by p = r/b(t). Finally, if we arrange the scale parameter 
b(t) according to 

the equation of motion for ip(r, p) becomes independent of b(t) 

i^(T, p) = (-± ^ + f ulp 2 + g tfY-W- 1 ) 4>(t, P) , (34) 
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i.e., we obtain a perfect scaling solution of the full field operator exactly in those cases 
where the effective metric is flat. Note that the implications of this property of the full 
field operator go far beyond the scaling of the hydrodynamic solution - we also obtain a 
perfect scaling of the quantum fluctuations to arbitrary order and for large wavenumbers 
(provided that the s-wave approximation is still valid, of course), where hydrodynamic 
(Thomas-Fermi) solution breaks down. It is not even necessary to assume the mean-field 
expansion, i.e., the scaling also applies to non-condensed bosons. 

Interestingly, all three examples (from D = 1 to D = 3) are potentially relevant for 
real physical systems: 

• D — 1 and N = 3: In quasi- ID Bose-Einstein condensates, one may obtain an 
effective coupling if the perpendicular size a z is comparable to the s-wave 
scattering length a s , see, e.g., (TUJ [TTJ [T2] . 

• D = 2 and N = 2: This is the usual quasi-2D Bose-Einstein condensate with 
quartic coupling, e.g., 

• D = 3 and N = 5/3: Even though the exponent N = 5/3 may seem somewhat 
unnatural, it appears in an effective description of a weakly interacting two- 
component Fermi gas in the BCS state, where the equation of state is determined 
by the Fermi energy E-p oc g 5 ^ 3 , and which therefore also shows scaling behaviour, 
see, e.g., |Hj. 

5. Horizon Analogues 

Apart from the question of whether it can be cast into a flat space-time form by means of 
an appropriate coordinate transformation or not, the emergence of a non-trivial effective 
metric in Eq. (fiT)|) suggests the application of concepts known from general relativity, 
such as horizons jU |2U 122] • Generally speaking, horizons correspond to a loss of causal 
connectivity, i.e., events beyond a horizon have no influence or cannot be influenced. 
The most prominent example is the event horizon of a (classical) black hole, beyond 
which everything is trapped forever (i.e., nothing can ever escape to infinity). 

On the other hand, in cosmology, two slightly different horizon concepts play a more 
important role due to the large-scale homogeneity and isotropy - the particle horizon 
and the apparent horizon. The particle horizon always refers to a chosen trajectory and 
indicates the border to the space-time region which cannot be reached by any signal 
starting from this trajectory or from where no signal can reach this trajectory. The 
apparent horizon depends on the chosen coordinates and is (roughly speaking) defined 
as the border beyond which all closed two-surfaces can either only expand or only 
contract w.r.t. the chosen coordinates. 

As it turns out, the concepts of the particle and the apparent horizon can be 
applied to expanding Bose-Einstein condensates. Let us start with the particle horizon 
and choose the trajectory r = — > p = for which the particle horizon can conveniently 
be determined using the metric in Eq. (|23|). The question is: which values of the co- 
moving coordinate p can a sound wave reach by starting at the origin p = at time 
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t and being described by a null line ds"l s = 0, i.e., how far can it travel? According to 
Eq. ([23)1 . this length can be calculated via the following integral (cf. also [3]) 

oo . oo . 

f I c N{t ) f dt 

A Phorizon (t) = J dt -^y = c N {t = 0) J &1+g(JV _ 1)/2(t/) • (35) 

If the above integral converges to a finite result Aphorizon(^), this corresponds to a particle 
horizon since no point beyond this co-moving coordinate can be reached by sound waves 
anymore. 

For a free expansion of the condensate V ext (t | oo) = 0, we obtain b(t f oo) oc t at 
late times according to Eq. (}3*3*|) and thus a horizon exists for all N > 1. The occurrence 
of a horizon might be puzzling if the metric in terms of the proper time r in Eq. (|24|) is 
flat - but this puzzle can be resolved by the observation that the proper time r reaches 
a finite value r(t f oo) < oo in the limit of arbitrarily late laboratory times in that case, 
cf. Eq. 

In contrast to the particle horizon, which necessitates the knowledge of the full 
future (or past), the apparent horizon can be identified by means of the configuration 
at a certain instant of time only. Instead of a trajectory as for the particle horizon, 
we have to choose a specific coordinate system (time slices) in order to define the 
apparent horizon. The coordinates (r, p) leading to a flat metric, for example, do of 
course not allow the introduction of an apparent horizon. On the other hand, since we 
observe the expanding condensate using the laboratory coordinate system (t, r), we shall 
choose these coordinates instead. For a spherically symmetric metric as in Eq. (|19j). the 
apparent horizon is determined by <?qo = 0' i- e -' where the velocity of the condensate v 
exceeds the speed of sound 

2 2 fcl-£>(JV-l)/2 

v o(t, ^horizon) = c jvW ¥ ^horizon (t) OC -. . (36) 



The very intuitive interpretation is that no sound wave can enter the region r < rh or izon 
from the outside. 

For perfect scaling N = 1 + 2/D and a freely expanding condensate, the apparent 
horizon settles down at late times (when b becomes constant) to a finite value 
^horizon (t T 00 ) = r horizon > 0) otherwise it may increase or decrease forever (depending 
on D and N). Note that both, the apparent and the particle horizon are never at rest 
w.r.t. the co-moving coordinate p but always decreasing. As a result, the wavelength 
of all sound modes exp{in ■ p} exceeds the horizon size at some point of time (horizon 
crossing) in view of the permanent stretching of modes due to the expansion of the 
condensate. After the wavelength exceeds the horizon size, the regions of higher and 
lower pressure cannot interact anymore and hence the modes stop oscillating (freezing of 
modes). As a result, comparing each mode with a harmonic oscillator, the momentum 
and its variance decrease drastically - and according to the Heisenberg uncertainty 
relation, the complementary variance must increase in order to compensate this. Ergo, 
the initial ground/vacuum state gets squeezed, which amplifies the quantum fluctuation 
in a certain direction. This admittedly rather intuitive picture (horizon crossing, 
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freezing, and squeezing) applies to both, the expanding universe (provided a horizon 
exists, such as during cosmic inflation) as well as expansing Bose-Einstein condensates 
in a very similar way and must of course be further supported by explicit calculations, 
see the next Section. 

Note that analogue horizons (and perhaps horizons in real gravity as well) are 
only low-energy effective concepts since, for very small wavelengths (below the healing 
length), the quantum pressure term becomes important leading to group and phase 
velocities exceeding the usual speed of sound. Furthermore, as the healing length 
increases during the expansion, the late-time limit in the derivation of the particle 
horizon is not strictly valid. Nevertheless, the main effects such as the amplification of 
quantum fluctuations are not modified drastically, see the next Section. 



6. Density-Density Correlations in quasi-2D 

According to the results of Section HJ the time-evolution of a quasi-two-dimensional 
condensate with the usual quartic coupling (i.e., N = D = 2) is particularly simple 
to solve due to its perfect scaling behaviour - which enables us to derive its dynamics 
exactly (e.g., including the quantum-pressure corrections). If we assume that the initial 
trapping potential is harmonic, the free expansion after switching off the longitudinal 
trap V\? xt (t > 0) = is fully described by the time-evolution of the scaling parameter 
according to Eq. (J33J) 



6(f) = VI + # • (37) 
The associated proper time r can be calculated from Eq. ()32|) 

, „ s arctana>nf , N 

r (t > 0) = . (38) 

u 

Since r quickly approaches a finite value at late times f f oo, the fluctuations freeze 
in and assume the value at proper time r = ti/{2ujq). In view of the perfect scaling 
discussed in Sec. EJ both the background and the phonon modes - expressed in terms of 
p and r - behave in the same way before and during the expansion, leaving the relative 
correlation function 

= (5Q(t,p)5g(t,f/)) = (5g(p)6g(p')) 
{P ' P) - g (t,p)g (t }P ') ~ q (p) 6 o(p') 1 ) 

unchanged. Consequently, the spectrum of the relative density contrast in co-moving 
coordinates after horizon crossing is given by the initial correlation spectrum. Besides 
the quasi-2D condensate with usual quartic coupling (N = D = 2), this also holds for 
all other cases where perfect scaling occurs. 

As discussed in Section the apparent horizon depends on the chosen time-slicing, 
i.e., the set of coordinates. (The particle horizon is independent of the set of coordinates, 
but requires the knowledge of the whole future evolution.) Co-moving coordinates, albeit 
advantageous to describe the evolution of the modes, are not suitable for describing the 
appearance of an apparent horizon. On the other hand, all observations are performed 
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in laboratory time t and coordinates r with the metric (|19|) showing the formation of 
an apparent horizon according to Eq. ([36)1 . For iV = 1 + 2/D, it follows 



rhorizon(t) = ~s = C iv(0)- — , (40) 



cn(o) , m yi + ^ 2 

= cjv(U)- ^— 

in the homogeneous case. While the trapping potential is still turned on (t < 0), the 
(apparent) horizon size rhorizon is infinite. Clearly, without any expansion, all points 
are causally connected via sound waves. After releasing the condensate, however, 
the horizon settles down very quickly {uj^t > 1) at a finite position (in laboratory 
coordinates) cn (0)/uo, given by the initial speed of sound cjv(O) and trapping frequency 

LO . 

Since the (apparent) horizon approaches a finte position in laboratory coordinates 
(t, r) and every phonon mode with a given wavelength A expands with the condensate 
cloud, each phonon mode will cross the horizon eventually. Equivalently, in terms 
of the co-moving coordinates, the horizon size decreases for all times (as long as the 
condensate is expanding) and the (co-moving) wavenumber n remains constant. As 
mentioned before, the initial correlation function of a quasi-2D condensate translates 
directly into the frozen density contrast after horizon crossing. 

In order to derive the density- density correlation function quantitatively, we 
quantise the phonon modes within the initial condensate Vq = using the approximation 
(centre of the BEC cloud) of a constant background density Qq = const. After a normal 
mode expansion, the Lagrange function (fT^jl reads 

- 1 _ 1.2 



4 V 



(41) 



where Vq denotes the quantisation volume. Note that, in contrast to Eq. (fTlj) . we did not 
neglect the quantum pressure term. In a straightforward manner, the phase and density 
fluctuations can be quantised via the introduction of creation (a h ) and annihilation 
operators (fife) which diagonalise the Hamiltonian. The quantised phase and density 
fluctuations read 



1 / k 2 



\ V Q uj k \ 4mpo 



TT \92D + ~r— J (ofc - 4) > (43) 



5S k 



\ Vq \ AmQo 
with the well-known Bogoliubov dispersion relation 

= V — + , 44 

where \x = (?2d£?o D is the (initial) chemical potential. 

For a perfectly scaling condensate, the two-point function in Eq. (p?Hj) yields the 
density-density correlations at all times, and in particular also the frozen correlations 
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(the density contrast) after horizon crossing. We can thus calculate the Fourier 
components 



in order to obtain the spectrum. 

For large wavelengths, the spectrum is linear in the (co-moving) momentum 
k in complete agreement to the curved space-time analogy. For large (co-moving) 
momenta k,, where the dispersion relation changes from linear to quadratic and 
the curved space-time analogy breaks down, the spectrum of the two-point correlation 
function approaches a constant value corresponding to a Dirac (^-contribution C(p, p') oc 
5(p— p') + . . .. Interestingly, after subtracting this local term, we obtain a scale-invariant 
spectrum C(k) oc 1/k 2 for large n - in analogy to cosmic inflation, where the spectrum 
is scale-invariant, too (in 3+1 dimensions it is 1/k 3 ). However, one should bear in mind 
that this scale-invariance occurs in a region where the curved space-time analogy breaks 
down. 

For addressing the question of whether this effect (amplification of the initial 
quantum fluctuations of the phon modes) can be measured, we estimate the order of 
magnitude of the relative density- density correlations C(k) for p ^ p' . Since the two- 
dimensional coupling constant, (72D; is proportional to the ratio of the scattering length, 
a s , and the perpendicular extension of the condensate, a z , we have typically 



Note, however, that this can only be seen as a rough estimate of the order of magnitude, 
since the two-point function is logarithmically divergent at small distances (where 
eventually the employed approximations break down). Remembering the conditions (0) 
in Section |21 we see that the relative density contrast must be significantly smaller than 
one - which was to be expected. However, the fluctuations could be on the percent-level 
(see Section |HJ) and thus might well be measurable. 

7. Density-Density Correlations in 3D 

For an expanding 3D condensate with quartic coupling the situation is more complicated, 
since it does not show the perfect scaling behaviour as in two dimensions. Therefore, it 
is not possible to infer the final (frozen) density (or phase) fluctuations from the initial 
state exactly as in the previous case. Fortunately, for long wavelengths, we may exploit 
the effective space-time analogy and derive the evolution using the tools known from 
cosmology. To this end, we transform onto co-moving coordinates and thus diagonalise 
the effective PGL-metric ()19|) . Note that, in three (spatial) dimensions, the background 
density scales as pit) = p(0)/6 3 (t). Assuming the background density to be smooth 
(Thomas- Fermi limit), the evolution equation for the scale factor - after switching off 




(45) 




(46) 
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the trap - reads, cf. 

1 ® = WY (47) 

Similar to the previous case, the scale factor accelerates (the interaction energy is 
transformed into kinetic energy) after switching off the trap, but quickly assumes a 
linear behaviour, b oc t (for t ^> \/ojq). 

The following calculations are most conveniently performed using the laboratory 
time but co-moving spatial coordinates. One advantage of the analogy to curved space- 
times is the inherited freedom of choosing suitable coordinates, where the independence 
of the final result follows automatically from covariance. The equation of motion of the 
field modes <p K with co- moving wavenumber k assumes the form 

[w + hm + ir)^ = «- («) 

Initially, for the trapped condensate, the damping term 3<p K b/b vanishes and the 
modes oscillate freely. During the expansion, the third term decreases whilst the 
second one increases and eventually dominates the latter one - in complete analogy 
to an over-damped oscillator, the modes freeze in. Because most phonon modes have 
frequencies much larger than the trapping frequency (e.g., [23-), this horizon crossing 
and freezing process happens during the period of linear expansion b(t) = at (with 
a ~ 0.82a;o). In view of the adiabatic theorem, the initial vacuum state for these 
(initially) rapidly oscillating modes is preserved (adiabatic vacuum) during the initial 
non-linear expansion. Thus the equation of motion describing the freezing in of modes 
can be simplified considerably by inserting b(t) = at 

(£ + 4i + 

and solved analytically in terms of Bessel functions [21] 

6 _ c o)l H m (l^l t -3/A i cW~h {2) (- r 3 / 2 ) (so) 

^"^^2/3^5/2* ^2/3^5/2* ) > ^ 

where Hyf denote the Hankel functions with the index v = 2/3 and are the 

corresponding integration constants. As one would expect, for early times, the modes 
<p K oscillate 0^' 2 ^ = izcUad^' 2 -*. It turns out that the Hankel functions have 
the proper asymptotic behaviour for early times such that the integration constants 
C^ 1 ' 2 -* can be replaced by (time-independent) creation and annihilation operators (in 
analogy to the previous Section) associated to the initial adiabatic vacuum state for the 
quantised phonon modes (fi K in the Heisenberg picture 

V 6 a 6 Vq 

where the pre-factor can be derived by imposing the canonical commutation relations 
for the conjugate variables <fi and 5q = <p/g obtained from the Lagrangian (|18|) for co- 
moving coordinates. Note that the argument of the Hankel functions coincides (up to a 
constant) with the proper time, r oc |a~ 5 / 2 t _3//2 , during the phase of linear expansion. 
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For late times t f oo, the modes <p n approach a constant value (due to horizon 
crossing and freezing) which can be derived inserting the asymptotic behaviour of the 
Hankel functions H^k [21] • The frozen late-time t j oo expectation value for phase- 
phase correlations reads 

r n ( t 1 / 3 S 4 / 3 

(<&) = J d 3 p cos( K ■ p)(0(p)0(O)) = |- [r (2/3)] 2 -^y^" 4 / 3 . (52) 

Employing Eq. we can calculate the density fluctuations, 5g = —<fi/g (in co-moving 
coordinates, where v vanishes), and obtain the spectrum of the relative density-density 
correlation function at late times 

where everything on the right-hand side (healing length etc.) refers to the initial state. 
Interestingly, one obtains the same spectra (A; -4 / 3 for the phase and /c +4//3 for the density 
fluctuations) for a condensate at rest after sweeping through the phase transition at 
g = by means of a time-dependent coupling g(t) [25] . 

In contrast to the quasi-2D case with perfect scaling, the above results rely on 
the curved space-time analogy and thus are only valid for wavelengths far above 
the healing length. Furthermore, in three dimensions, the healing length scales like 

= £(0)& 3 / 2 (t) and grows faster (again in contrast to the quasi-2D case) than the 
wavelengths X(t) = X(0)b(t) in laboratory coordinates. In order to obtain a rough 
estimate of the maximum amplification effect, we insert the maximum co- moving 
wavenumber k = 1/ 1 ^{a/uj^) 1 / 4 of the phonon modes whose wavelength exceeds the 
healing length during the relevant period of their evolutions (i.e., until horizon-crossing 
and freezing) 

, c (K) 4A[r(i/3)] W / 4 r—^r 3/4 



cu 5 \- 3 / 4 



« 30.3^0 -i . (54) 

Here = co/^(0) denotes the (initial) frequency of the phonons with the (co- moving) 
wavenumber l/£(0) and a = au>o where a ~ 0.82 is a dimensionless constant. In 
contrast to the quasi-two-dimensional case, the size of the fluctuations now depends on 
the diluteness parameter Ja^Qo (instead of the ratio a s /a z ), which must be small for the 



employed approximations to apply ya^g <^ 1. Furthermore, the ratio uj^/ujq = Cqk/ujq 
should be large uj k /ujq ^> 1 if we want a sufficiently long period of linear expansion 
b(t) = at. However, the smallness of these parameters can partly be compensated by 
the numerical pre-factor 30.3 such that the final effect can be on the percent-level (see 
the next Section). 

The spectrum and size of the frozen-in correlations C3D could be measured by 
obtaining a density map of a slice of the condensate (tomographic imaging). By making 
a projective (e.g., absorption) image instead, one would average over the fluctuations 
in the (spatial) direction of projection. While the spectrum remains unchanged, this 
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averaging yields an additional suppression by a factor of order £/l z with l z denoting the 
transversal (i.e., in direction of projection) extension of the condensate. 

8. Summary 

The objective was to investigate and to exploit the analogy between phonons in 
expanding Bose-Einstein condensates on the one hand and quantum fields in an 
expanding universe on the other hand. The advantages of this curved space-time 
analogy are two-fold: Firstly, it facilitates the application of all the concepts, tools, 
and effects known from general relativity (such as horizons) and thereby fosters a 
better understanding of condensed-matter phenomena in terms of universal geometrical 
concepts. Secondly, this analogy enables a theoretical and possibly experimental 
investigation of exotic quantum effects known from cosmology. 

After a discussion of the effectively lower-dimensional behaviour of strongly 
constrained BECs in Section |21 the emergence of an metric governing the propagation of 
low-energy phonons is discussed in Section EJfor arbitrary power-law self-coupling in any 
dimension. In complete analogy to cosmology, the introduction of co-moving coordinates 
is advantageous for the description of many phenomena. It turns out that the full 
field operator (in a harmonic trapping potential) possesses perfectly scaling solutions 
exactly if the effective metric is flat in terms of the co-moving coordinates. In a freely 
expanding Bose-Einstein condensate, an effective sonic horizon is formed, i.e., two points 
at fixed (co-moving) spatial positions whose distance exceeds the horizon size cannot 
be connected anymore by (low-energy) phonons, cf. Section 03 The formation of this 
effective horizon implies the amplification of the quantum fluctuations (horizon-crossing 
and freezing) already known from cosmology, which has been calculated explicitely for 
quasi-2D and 3D condensates with quartic coupling in Sections El and [3 

In order to obtain an explicite estimate for the size of the derived effect, let us 
consider a condensate of 10 5 sodium atoms inside a highly anisotropic trap with the 
trapping frequencies uj±/2tt = 790 Hz and uj\\/2it = 10 Hz, cf. [26\. The thickness of 
the disk, a z = 0.746 /im, is smaller than the healing length, £ = 1.34 /xm, but still 
much larger than the scattering length, a s = 2.8 nm, complying with the hierarchy of 
scales derived in Sec. El To estimate the order of magnitude of the density-density 
correlations, we consider a mode with (co-moving) wavenumber k = 27r/£. For the 
relative density contrast inside a volume of size £ 2 we obtain C(k)/^ 2 = 1.79%. For 
the modes in the linear regime of the spectrum, horizon crossing and thus freezing 
in occurs shortly after the trap is switched off. Initially the horizon is at infinity, 
but very quickly settles at rh or izon(£ T °°) = 3.28 /im. Hence already after one e-fold, 
when the radius of the atom cloud has increased by a factor e ~ 2.7, all modes with 
wavelength larger than £ are frozen in and the fluctuations are transformed into a density 
contrast. With the correlation function depending on the ratio a s /a z , the effect can 
be enhanced by confining the condensate more tightly in the perpendicular direction. 
However, if a s /a z is not small, the condensate is no longer described by the usual Gross- 
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Pitaevskii Lagrangian and our analysis does not apply anymore in this form (e.g., the 
perfect scaling breaks down). For such a case, the (initial) correlation spectrum (with 
appropriate non-quartic coupling, cf. [121 EZ|) was calculated in Ref. |2Z| and differs 
from our result in Eq. (j45jl . 

As an example for the 3D case, let us consider 87 Rb atoms inside a spherically 
symmetric trap. For a condensate consisting of 10 7 atoms and a trapping frequency 
ujq/2tt = 200 Hz, which are potentially realisable parameters, cf. [2*S] . the Thomas-Fermi 
radius is R = 12.2 /im. The minimal frequency of phonons uWi = 2ttc s /R = 5582 Hz is 
well above the trapping frequency uq. Hence the freezing of modes takes place during 
the period of linear expansion and the approximation b(t) = at in the evolution equation 
of the phonon modes is justified. A rough estimate of the maximum effect according to 
Eq. flUJ) yields 



i.e., a potentially measurable effect 

So far, we have assumed zero temperature. Of course, initial thermal fluctuations 
would - provided that there is no thermalisation during the expansion of the condensate 
- also be amplified by basically the same mechanism as the quantum fluctuations. In 
order to ensure that the quantum fluctuations are larger than the (initial) thermal 
fluctuations, the temperature must be small enough such that the thermal occupation 
of the phonon modes under consideration is negligible. For the parameters used above, 
the relevant temperature scales are T(l/£) = 0.1 nK for the Q2D and T(k) = 3.9 nK 
for the 3D condensate, respectively. By increasing the particle number and/or the 
trapping frequency, the requirements regarding the experimental temperature can be 
relieved. (Another way to discriminate between thermal and quantum fluctuations is 
the spectrum.) 

In summary, expanding Bose-Einstein condensates facilitate the experimental 
simulation of exotic effects of quantum fields in curved space-times. On the other 
hand, the amplification mechanism under consideration sets the ultimate quantum limit 
of accuracy for time-of-flight experiments: No matter how smooth the initial cloud 
can be prepared, the frozen and amplified quantum ground-state fluctuations generate 
noticeable density perturbations. 
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